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I. INTRODUCTION 

One of the most striking aspects of quantum coherence 
in a quantum many-body system is the occurrence of en- 
tanglement, namely the realization of a superposition of 
many-body states that cannot be factorized into a prod- 
uct of single-particle wave functions. An entangled state 
possesses correlations that cannot be accounted for by 
classical-like quantities; for instance an entangled state 
might not show any form of classical order, and nonethe- 
less be at the same time strongly correlated. The pos- 
sibility of a local description of such state is partially 
or completely lost, depending on the degree of entangle- 
ment contained in the state. In particular, the non-local 
nature of special collective quantum states is the funda- 
mental ingredient that allows quantum communication 
protocols and quantum computation algorithms^ to out- 
perform their classical counterparts. 

The idea of entanglement as a resource naturally 
demands a systematic investigation of which quantum 
many-body systems are able to display sizable entangle- 
ment in a controllable way. An intriguing perspective 
is that pure-quantum correlations are strongly enhanced 
when a system undergoes a quantum phase transition 
(QPT)2i analogously to what classical correlations do at 
a thermal phase transition. Indeed quantum fluctuations 
show up at all length scales at a quantum critical point. 
In what sense quantum correlations, and thus entangle- 
ment, 'diverge' at a QPT has been the subject of investi- 
gation in several recent studiesi^*^, although the result- 
ing picture is still controversial. In this work we consider 
a quite general S = \/2 quantum spin system displaying 
a field-induced QPT in arbitrary dimensions. In particu- 
lar, we discuss how the behavior of entanglement at and 
around the quantum critical point shows strong features 
providing new insight in the drastic change of the sys- 
tem's ground state under the effect of strong quantum 
fluctuations. 



II. THE MODEL 

A very general example of a quantum phase transition 
in spin models is offered by the the antiferromagnetic 



XYZ model in a field. The model Hamiltonian reads: 
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where J > is the exchange coupling, the sum (ij) runs 
over the nearest neighboring sites of a bipartite lattice 
with coordination number z, and h = gfi-QH/J is the re- 
duced magnetic field. In the following we perform the 
canonical transformation S^'^ {~iyS^'^ on Eq. 
so that the relevant correlations along the x and y axes 
are ferromagnetic. The parameters < Aj,, < 1 con- 
trol the anisotropy of the system. In the most general 
case of Aj, ^ 1 the Zeeman term in Eq. does not com- 
mute with the rest of the hamiltonian. This property 
is at the core of the field-driven quantum phase transi- 
tion occurring at a critical field /ic(A), which separates 
a Neel-ordered phase {h < he) from a partially-polarized 
disordered phase {h > he). When h < he the field favors 
long-range magnetic order along the x-axis. This order 
disappears at the critical field he, where magnetic cor- 
relations along X become short-ranged, while quantum 
fluctuations prevent the spins from being fully polarized 
along the field2i2iii. 

The case A^ = reproduces the XY model in a trans- 
verse field, which is exactly solvable in one dimension^, 
and whose entanglement properties have been the sub- 
ject of several recent investigation3MiS4&. When A^ ^ 0, 
and/or in higher dimensions, the model is no longer 
exactly solvable. The one-dimensional case has been 
indeed investigated within approximate analytical and 
numerical approachesSiifliii. A renewed interest in the 
model stems from the experimental in-held studies on 
the quantum spin chain compound Cs2CoCl4^^, display- 
ing a strong planar anisotropy, Aj, « 0.25, A^ w 1, and 
J w 0.23 meV. 



III. LINEAR CHAIN: CRITICAL POINT AND 
FACTORIZED STATE 

Motivated by the existing theoretical and experimental 
results, we first concentrate on the case of a linear chain, 
and on the parameter range < Aj, < 1, A^ = 1, which 
defines the XYX model in a field^^. The qualitative be- 
havior of the more general XYZ model is expected to be 
close to that of the XYX model in a field, since the two 
models share the same symmetries. 
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FIG. 1: Upper panel. Ground state phase diagram of the 
XYX model in a field. Mean-field (MFA) results are taken 
from Ref-— , the factorizing field h{ from Ref>^. Lower panel. 
Quantum critical behavior of x— and z— magnetizations and 
correlation length for the model with Aj, = 0.25, L = 100, 
13 = 200. The factorizing field is indicated by a dashed line. 
The arrow indicates the quantum critical point (QCP). 

The analysis of the theoretical model is performed 
via Stochastic Series Expansion (SSE) Quantum Monte 
Carlo (QMC) simulations^^, based on a modified version 
of the directed-loop algorithm to account for the low sym- 
metry of the Hamiltonian^^. Chains of various lengths, 
L = 40, 120, have been considered, at an inverse tem- 
perature (3 — 2L high enough to mimic the T — be- 
haviour for each lattice size. 

The left panel of Fig. shows the ground-state phase 
diagram of the one-dimensional XYX model in the Ay — 
h plane. The quantum Monte Carlo data confirm that 
the transition belongs to the universality class of the ID 
transverse- field Ising model (or of the 2D Ising model)^, 
and are in very good agreement with predictions from a 
mean-field treatment of the Hamiltonian^. On the right 
panel of Fig. ^ we can clearly see the signatures of the 
transition in the critical behavior of the magnetizations 
along X (estimated through the asymptotic value of the 
spin-spin correlator as = \{Sf S^_^_^^2)\^^^) along 
z, and the divergence of the correlation length along x. 



Remarkably, none of these standard magnetic observ- 
ables bears signatures of the second striking feature of 
the model, namely the occurrence of an exactly factor- 
ized state for a field /if (Ay) lower than the critical field 
he- The factorizing field reads h{ — •\/2(l + Ay) in the 
case of the XYX model. The factorized state has the 
form I*) = <S>f=i IV'j) with IV-j) = cos6l| t) +sin6'e*'^^- 1 i) 
and (jjj = (1 + (-l)J>/2, e = cos-i ^{l + Ay)/2. This 
corresponds to a configuration in which the spins have 
perfectly staggered x components but also cant out of 
the xy plane by an angle 9. The occurrence of such a 
factorized ground state is particularly surprising if one 
considers that we are dealing with the 5 = 1/2 case, 
characterized by the most pronounced effects of quan- 
tum fluctuations. However, in the class of models here 
considered, such fluctuations are fully uncorrelated^i at 
h — hi, thus leading to a classical- like ground state. 

IV. ENTANGLEMENT ESTIMATORS 

The study of entanglement properties in the XYX 
model turns out to be very insightful, and it has the 
unique feature of unambiguously detecting both the fac- 
torized state and the quantum critical point. Our study 
has focussed on the entanglement of formatior}^ which 
is quantified through the one-tangle and the concur- 
rence. The one-tangle^^J^ is an estimate of the T = 
entanglement between a single site and the remainder 
of the system. It is defined as ti = 4detp(^\ where 
= (/ -f Af"(T")/2 is the one-site reduced den- 
sity matrix, M" — (S""), are the Pauli matrices, and 
a = x,y, z. In terms of the spin expectation values Af 
Ti takes the simple form: 

n = 1-4^(M")2. (2) 

The concurrence^- quantifies instead the pairwise entan- 
glement between two spins at sites i, j both at zero and 
finite temperature. For the model of interest, in absence 
of spontaneous symmetry breaking (M^ = 0) the con- 
currence takes the form.^-' 

=2max{0,C«,Cf } , (3) 

where 

- 9tf-l + W-9m, (4) 

Cg^ = W +9ff\- ^(1+9!!^ -{M^y , (5) 
with = 

V. RESULTS 

The QMC results for the model Eq. (P) with = 0.25 
are shown in Fig. [3 where we plot ri , the sum of squared 
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FIG. 2: Upper panel. One-tangle ri and sum of squared con- 
currences T2 as a function of the applied field for the S = 1/2 
XYX model with Ay = 0.25, L = 100 and /3 = 200. Inset: 
contributions to the concurrence between j-th neighbors; full 
symbols stand for Cf'^'^j , open symbols for C^'^^^j. The dashed 
line marks the critical field he- Middle and lower panel. Spin- 
spin correlator gf^+r compared with the concurrence Ci,i+r 
for various field values around the critical value he — 1.605(2). 
Other simulation parameters as in the upper panel. 
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and, in the inset, Ci_i+„ for n = 1,2,3. The following 
discussion, although directly referred to the results for 
Ay — 0.25, is actually quite general and applies to all 
the other studied values of Aj,. 

Unlike the standard magnetic observables plotted in 
Fig. ^ the entanglement estimators display a marked 
anomaly at the factorizing field, where they clearly van- 
ish as expected for a factorized state. In particular, the 
power of the entanglement estimators is to rigorously un- 
veiling the factorized nature of the ground state. It can 
indeed be easily shown that the ground state is factorized 
if and only if the one-tangle ti vanishes for each spin. 
When the field is increased above /if, the ground-state 



FIG. 3: Entanglement ratio T2/T1 as a function of the field 
for Ay — 0.25 and /3 — 2L. Inset: zoom on the critical region. 



entanglement has a very steep recovery, accompanied by 
the quantum phase transition at he > hf. The system re- 
alizes therefore an interesting entanglement switch effect 
controlled by the magnetic field. 

As for the concurrence terms Eqs. Q),©, the factor- 
izing field divides two field regions with different expres- 
sions for the concurrence: 



(1) 



for h < hf , 
for h > h[ , 



(7) 
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whereas C^^ = C^-^^ — at h = hi. In presence of 
spontaneous symmetry breaking occurring for h < he, 
the expression of the concurrence is generally expected 
to change with respect to Eqs. 10}, as extensively dis- 
cussed in Ref^ . For the model under investigation, the 
expression of the concurrence stays unchanged when the 

conditionSS cl^^ < cl^^ is satisfied, i.e. for h > hr. This 
•J , -J 

means that our estimated concurrence is accurate even 
in the ordered phase above the factorizing field; in the 
region < h < h{ it represents instead a lower bound 
to the actual T = concurrence. Alternatively it can be 
regarded as the concurrence for infinitesimally small but 
finite temperature. 

A naive analogy between classical and quantum cor- 
relations would lead to the expectation that the range 
of pairwise entanglement, expressed through the concur- 
rence, is critically enhanced at a quantum phase tran- 
sition, reflecting the divergence of the length scale for 
quantum effects. Indeed this is clearly not the case as 
shown in the right panel of Fig. [21 where the behavior of 
the spin-spin correlator g^^ as a function of the distance 
is contrasted with that of the concurrence. We observe 
that, while the correlator becomes long-ranged below the 
critical field and even completely flat at the factorizing 
field, the concurrence remains short-ranged when pass- 
ing through the transition, and it basically vanishes after 
four lattice spacings. 



To clarify this issue, one has to consider the special na- 
ture of quantum correlations. At variance with classical 
correlations, entanglement has basically to satisfy a sum 
rule, or, as it is often phrased, is subject to a constraint 
of monogamy. This means that the more partners a spin 
is entangled with, the less entanglement will be shared 
with each partner. Moreover, entanglement is present not 
only in the form of pairwise correlations as in the case 
of the concurrence, but also in the form of n-spin corre- 
lations with n > 2, which, unlike classical correlations, 
can be completely independent of pairwise correlations. 
This is exemplified for instance by the maximally en- 
tangled Greenbergcr-Horne-Zeilinger (GHZ) state for n 
spins |GHZ) = (I tt - T) + Ui - i))/ A on which the 
concurrence between any two spins is vanishingl- . There- 
fore, according to monogamy, pairwise entanglement and 
n-wise entanglement with n > 2 are mutually exclusive, 
a condition which is absurd from the point of view of 
classical correlations. 

The mathematical expression of the monogamy con- 
straint for the entanglement is provided by the Coffman- 
Kundu-Wootters (CKW) conjecture^^, which indeed rep- 
resents an approximate sum rule for entanglement corre- 
lations. Such conjecture states that T2 < ti, and, al- 
though it can be rigorously proved only in the case of 
three spins, it has always been verified so far in the case 
of an arbitrary number of spins^^. Indeed Fig. |21 shows 
that the CKW conjecture is verified also in the case of 
our model. 

In particular, we interpret the entanglement ratio R — 
T2 /ti as a measure of the fraction of the total entangle- 
ment stored in pairwise correlations. This ratio is plotted 
as a function of the field in Fig. O As the field increases, 
we observe the general trend of pairwise entanglement 
saturating the whole entanglement content of the system. 
But a striking anomaly occurs at the quantum critical 
field he, where R displays a very narrow dip. According 
to our interpretation, this result shows that the weight 
of pairwise entanglement decreases dramatically at the 
quantum critical point in favour of multi-spin entangle- 
ment. Indeed, due to the monogamy constraint, multi- 
spin entanglement appears as the only possible quantum 
counterpart to long-range spin-spin correlations occur- 
ring at a quantum phase transition. The divergence of 
entanglement range has to be interpreted in the sense 
that, at a quantum phase transition, finite 71-spin entan- 
glement appears with n — > cx) at the expense of the pair- 
wise one. A unique estimator for the n-spin entanglement 
with n > 2 has not been found yet, so that justifying the 
above statement on a more quantitative level is problem- 
atic. Nonetheless strong indications of the relevance of 
n-spin entanglement with large n at a quantum critical 
point are given in Ref.^. Finally, the above results evi- 
dence the serious limitations of concurrence as a reliable 
estimate of entanglement at a quantum critical point. In 
turn, we propose the minimum of the entanglement ra- 
tio i? as a novel estimator of the quantum critical point, 
fully based on entanglement quantifiers. 
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FIG. 4: One-tangle and sum of squared concurrences as a 
function of the field for the XYX model on the two-leg ladder 
with L = 40 on each leg and Ay = 0. Inset: entanglement 
ratio and z-magnetization as a function of the field. The 
shaded area marks the quantum disordered region. All lines 
are guides to the eye. 



VI. TWO-LEG LADDER 

We extend the above analysis of entanglement to the 
less investigated case of the XYX model on a two- 
leg ladder, thanks to the fact that the SSE quantum 
Monte Carlo approach is easily generalizable to any bi- 
partite lattice. In this case the hamiltonian of Eq.(^) is 
still expected to display a field-driven quantum phase 
transition, whose existence is independent of the lat- 
tice geometry and it is instead uniquely due to the non- 
commutativity between the exchange term and the Zee- 
man term. There is instead no specific reason to expect 
that the model still displays a perfectly factorized state 
for geometries other than that of the linear chain. 

Fig. 21 shows the one-tangle and the sum of squared 
concurrences for the XYX model on the two-leg ladder 
in the case — 0. We notice that, for such a strong 
anisotropy, the two-leg ladder does not display a gapped 
Haldane phase, as shown by the finite value of the mag- 
netization for low fields (inset of Fig. 01 . Remarkably, 
the qualitative behavior of both ri and T2 as a function of 
the field is the same as in the case of a single chain, and 
in particular the vanishing of the one-tangle signals the 
rigorous existence of a factorized state for a lower field 
than the critical one. The CKW conjecture is verified for 
all the field values considered, and therefore the entan- 
glement ratio T2 /ti can be still interpreted as a measure 
of the fraction of entanglement stored in pairwise cor- 
relations. This quantity, plotted in the inset of Fig. 01 
again displays a deep minimum corresponding to the in- 
flection point of the uniform magnetization and therefore 
marking the quantum critical point. 
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VII. CONCLUSIONS 

In this paper we have shown that entanglement esti- 
mators give a precious and novel insight in the ground 
state properties of lattice S = 1/2 spin systems. In the 
case of anisotropic spin chains with a field-driven quan- 
tum phase transition, we have shown that the quantum 
critical point can be detected through a narrow dip in the 
pairwise-to-global entanglement ratio. Moreover, unlike 
the more conventional magnetic observables, entangle- 
ment estimators are able to single out the occurrence of 
an exactly factorized state in these systems. The use 
of quantum Monte Carlo techniques naturally allows to 
extend this analysis to different lattice geometries. In 
particular we have discussed here the case of a two-leg 
ladder, in which the calculation of entanglement estima- 
tors remarkably shows the existence of a factorized state 
below the quantum critical point. The entanglement ra- 
tio displays again a minimum at the critical field, con- 
firming the generality of this feature of the entanglement 



behavior as a signature of a quantum phase transition. 
Finally, the proximity of a quantum critical point to the 
factorized state of the system gives rise to an interest- 
ing field-driven entanglement-switch effect. This suggests 
that many-body effects, driven by a macroscopic param- 
eter as an applied field, are a powerful tool for the control 
of the microscopic entanglement in a multi-qubit system. 
The application of this kind of concepts in the design of 
quantum computing devices looks therefore appealing. 
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